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ABSTRACT 

Multiple-input multiple-output (MIMO) radar systems have been 
shown to achieve superior resolution as compared to traditional radar 
systems with the same number of transmit and receive antennas. 
This paper considers a distributed MIMO radar scenario, in which 
each transmit element is a node in a wireless network, and investi- 
gates the use of compressive sampling for direction-of-arrival (DOA) 
estimation. According to the theory of compressive sampling, a sig- 
nal that is sparse in some domain can be recovered based on far 
fewer samples than required by the Nyquist sampling theorem. The 
DOA of targets form a sparse vector in the angle space, and there- 
fore, compressive sampling can be applied for DOA estimation. The 
proposed approach achieves the superior resolution of MIMO radar 
with far fewer samples than other approaches. This is particularly 
useful in a distributed scenario, in which the results at each receive 
node need to be transmitted to a fusion center for further processing. 

Keywords: compressive sampling, compressive sensing, 
MIMO radar, DOA estimation 

1. INTRODUCTION 

Unlike a conventional transmit beamforming radar system that 
uses highly correlated waveforms, a multiple-input multiple- 
output (MIMO) radar system transmits multiple independent 
waveforms via its antennas |[TJ - (5J . A MIMO radar system is 
advantageous in two different scenarios J5J. In the first one, 
the transmit antennas are located far apart from each other 
relative to their distance to the target. The MIMO radar sys- 
tem transmits independent probing signals from decorrelated 
transmitters through different paths, and thus each waveform 
carries independent information about the target. Therefore, 
the MIMO radar system can reduce the target radar cross sec- 
tions (RCS) scintillations and provide spatial diversity. In the 
second scenario, a MIMO radar is equipped with M t transmit 
and M r receive antennas that are close to each other relative 
to the target, so that the RCS does not vary between the differ- 
ent paths. In this scenario, the phase differences induced by 
transmit and receive antennas can be exploited to form a long 
virtual array with M t M r elements. This enables the MIMO 
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radar system to achieve superior spatial resolution as com- 
pared to a traditional radar system. In this paper we consider 
the second scenario. 

Compressive sensing (CS) has received considerable at- 
tention recently, and has been applied successfully in diverse 
fields, e.g., image processing [6] and wireless communica- 
tions Q. The theory of CS states that a K -sparse signal x 
of length N can be recovered exactly from 0(K log 7Y) mea- 
surements with high probability via linear programming. Let 
denote the basis matrix that spans this sparse space, and 
$ a measurement matrix. The convex optimization problem 
arising from CS is formulated as follows: 

min ||s||i, subject to y = $x = t&Vl/s (1) 

where s is a sparse vector with K principal elements and the 
remaining elements can be ignored; $ is an M x N matrix 
incoherent with and M < N. 

The application of compressive sensing to a radar system 
was investigated in JH), J9) and iflOl . In J8), it was demon- 
strated that the CS method can eliminate the need for match 
filtering at the receiver and has the potential to reduce the re- 
quired sampling rate. In the context of Ground Penetrating 
Radar (GPR), J9) presented a CS data acquisition and imaging 
algorithm that by exploiting the sparsity of targets in the spa- 
tial space can generate sharper target space images with much 
less CS measurements than the standard backprojection meth- 
ods. Also the sparsity of targets in the time-frequency plane 
was exploited for radar in 1101 . In the context of communica- 
tion, |[TT1 proposed the direction of arrival estimation (DOA) 
estimation using CS. In ifTTI . the basis matrix "J/ is formed by 
the discretization of the angle space. Since the signal sources 
were assumed to be unknown, the basis matrix was approxi- 
mated based on the signal received by a reference vector. That 
signal would have to be sampled at a very high rate in order 
to construct a good basis matrix. 

In this paper, we extend the idea of l8l- lfTT1 to the prob- 
lem of DOA estimation for MIMO radar. Since the number of 
targets is typically smaller than the number of snapshots that 
can be obtained, DOA estimation can be formulated as the re- 
covery of a sparse vector using CS. Unlike the scenario con- 
sidered in IfTTI . in MIMO radar the transmitted waveforms are 
known at each receive antennas, so that each receive antenna 
can construct the basis matrix locally, without the knowledge 



of the received signal at other antennas. Further, radar sys- 
tems often suffer from interference due to jammers. Jammer 
suppression is investigated here by exploiting the uncorrelat- 
edness of the transmitted waveforms with the jammer signal 
in order to design the measurement matrix. We provide an- 
alytical expressions for the signal-to-jammer ratio (SJR) for 
the proposed approach. We also provide simulation results 
to show that the proposed approach can accomplish super- 
resolution in MIMO radar systems by using far fewer sam- 
ples than existing methods, such as Capon, amplitude and 
phase estimation (APES) and generalized likelihood ratio test 
(GLRT) Q ■ This is very significant in a distributed scenario, 
in which the receive nodes would need to transmit the locally 
obtained information to a fusion center. For such systems, we 
show that the proposed approach can enable each node to ob- 
tain a good DOA estimate independently. Further, it requires 
much less information to be transmitted to a fusion center, 
thus enabling savings in terms of transmission energy. 

2. SIGNAL MODEL FOR MIMO RADAR 

We consider a MIMO radar system with Al t transmit anten- 
nas and M r receive antennas. For simplicity, we assume that 
targets and antennas all lie in the same plane. Let us denote 
the locations in rectangular coordinates of the i-th transmit 
and receive antenna by [x\, yf) and (x[, y\), respectively. We 
assume that all transmit and receive node locations relative 
to some reference point are known to each node in the net- 
work. In a clustered system this information pattern may be 
achieved via a beacon from the cluster-head (3J. The loca- 
tion of the fc-th target is denoted by (dk,0k), where dk is the 
distance between this target and the origin, and 9k is the az- 
imuthal angle, which is the unknown parameter to be esti- 
mated in this paper. 

Under the far-field assumption dk 3> \J {x\) 2 + (y*) 2 
and dk 3> J (x 7 -) 2 + (y[) 2 , the distance between the ith 
transmit/receive antenna and the fc-th target d\ k ld r ik can be 
approximated as d*£ r f» dk - i] 1 / 7 (6 k ), where ril /r (9 k ) = 
x\ /r cos(0 k ) + y\ /r sin(0 fe ). 

Let Xi in) denote the discrete-time waveform transmitted 
by the i-th transmit antenna. Assuming the transmitted wave- 
forms are narrowband and the propagation is non-dispersive, 
the received baseband signal at the fc-th target equals |4) 

Mi 
i=l 

(2) 

for fc = 1,..,,K, where x(n) = \x\(n), ...,XM t ( n )] T an d 
v(0fc) = [e i2 x^ ek \ e^riitiOk^T is the steer j n g vector. 

Due to reflection by the target, the Z-th antenna element 
receives 

K 

zt(n) = e- j ^"y k (n) + ei(n), l = l,...,M r (3) 
fc=i 



where e; (n) represents independent and identically distributed 
(i.i.d.) Gaussian noise with variance a 2 . 

On letting L denote the number of snapshots, and placing 
zi(n),n = 0, L — 1 in vector z; we have 

K 

zi=J2 (2<W (efc)) /?fcXv(# fc ) + e* (4) 
fe=i 

wherey fe = [j/fc(0), . . . , y k (L - 1)] T , e z = [ej(0), . . .,ei(L- 
l)] T andX= [x(0),...,x(i-l)] T . 

By discretizing the angle space as a = [ai, . . . , ajv], we 
can rewrite © as z; = J2 n =i e , ^ ,, ''' a "'stiXv(a n ) + e;, 
where N >> L, and 

f e^^ dk (3k if there is target at a n 
n \ otherwise 

3. COMPRESSIVE SENSING FOR MIMO RADAR 

Assuming that there exists a small number of targets, the DOAs 
are sparse in the angle space, i.e., s = [si, . . . , sjv] is a sparse 
vector. A non-zero element with index j in s indicates that 
there is a target at the angle oij . By CS theory, we can con- 
struct a basis matrix "f / for the Z-th antenna as 

= [e^"' (ai) Xv(ai), . . . , e^^ (Qjv) Xv(ajv)] . (5) 

Ignoring the noise, we have z; = ^fS. Then we measure 
linear projections of the received signal at the l-th antenna as 

ri = $/z ; = $/*/s (6) 

where $/ is an M x L random Gaussian matrix which has 
small correlation with Placing the output of N r receive 
antennas, i.e., ri, r]y r , in vector r we have 

r = 6s, 1 < Nr < M r . (7) 

and the structure of can be easily inferred based on ((6]). 
Therefore, we can recover s by applying the Dantzig selector 
to the convex problem in (f7]i as in IPT21 : 

s = mm||s||i s.t.||e H (r-0s)|| oo <^. (8) 

According to fPHl . we can recover the sparse vector s with 
very high probability to select /i = (1 + t~ 1 )^/2\ogNa 2 , 
where t is a positive scalar. 

4. PERFORMANCE ANALYSIS IN THE PRESENCE 
OF A JAMMER SIGNAL 

In this section, we analyze the effects of a jammer signal on 
the performance of DOA estimation for MIMO radar using 



CS in terms of the SJR. In the presence of a jammer signal, 
the received signal of the l-th receive antenna is given by 



fe=i 



-*ie -J '*< d-, ' r P ) >j8b + Siei 



(9) 



The location of the jammer is denoted by (d, 9), and (3 and b 
denote respectively the reflection amplitude and waveforms 
of this jammer. Since b is uncorrelated with the transmitted 
waveforms X, the effects of the jammer signal are similar to 
those of the addictive noise. Let A; = •J'f^;. Given the the 
transmit and receive node locations, the average power of the 
desirable signal P s (I) over the transmit waveforms is 

K 

P S {1) = E{Y, e^ 2 ^- d ^ - W (»* f3* k p k , 



k,k' = l 



pi(k,k') 



K 



: v g (fl fc )X H A ; Xv(0 fc Q } = E{J2 \Pk\ 2 Rkk} 



k=l 



Ci(0 



■E{J2pi(k,k')P* k (3 k ,R kk ,} 

k^k' 



(10) 



c 2 (i) 



where C±(l) is the dominant term and C*2(0 can be ignored 
when the number of transmit antennas is sufficiently large. 
Therefore, P s (l) can be approximated by C\(l). To sim- 
plify the derivation of the SJR, we assume that the transmitted 
waveforms are independently generated, orthogonal quadra- 
ture phase shift keyed (QPSK) sequences with X H X = Im^ 
b H b = 1 and &i$>f = 1m- Then P s (l) is approximated by 



K 

E 

k=l 



K 



^| 2 tr(A ; )M*/i=^El^| 2 

fc=i 



(11) 



Similarly, the average power of the jammer interference 
over the jammer waveforms is given by 



P 3 (l) 



E{ (e~^ ( d -rt(»))p) ( e -i ¥ (d-vl(e))py 
xb ff A z b} = \/3\ 2 M/L . (12) 



From these two expressions, the SJR becomes 



SJR 



(13) 



Since the jammer signal is uncorrelated with the transmit- 
ted signal, the SJR can be improved by correlating the jammer 
signal with the transmitted signal. Combining this with CS, 
the measurement matrix in (|6]l is modified as = $;X H . 
Moreover, since is a Gaussian random matrix, is still 



Gaussian; therefore it satisfies the restricted isometry prop- 
erty (RIP) and is incoherent with thus guaranteeing a sta- 
ble solution to ([8}. Based on $/, the SJR can be obtained as 



SJR 



(14) 



Generally, the SJR can be improved by a factor of Lj M t using 
$i since L ^> M t . ( TBI indicates that the increase in L will 
improve the DOA estimates. However, more calculations are 
required by £i-norm minimization due to the increase in the 
size of the basis matrix, and the time duration of the radar 
pulse needs to be longer as well. 

The proposed method is especially advantageous in a dis- 
tributed MIMO radar system in which the receive elements 
are randomly distributed. In particular, many fewer measure- 
ments are required to be sent to the base station or fusion 
center (FC) in this situation than are need by conventional 
methods. As simulation results show (see Section[5]), the pro- 
posed method can yield good performance even using a single 
receive antenna. With a good initial estimate of DOA, the re- 
ceive nodes can adaptively refine their estimates by construct- 
ing a higher resolution basis matrix "f/ around that DOA. Re- 
stricting the candidate angle space, may reduce the samples 
in the angle space that are required for constructing the basis 
matrix, thus reducing the complexity of the t\ minimization 
step. On the other hand, the resolution of target detection can 
be improved by taking the denser samples of the angle space 
around the intimal DOA estimate. 

5. SIMULATION RESULTS 

In this section, we consider a MIMO radar system with the 
transmit/receive antennas randomly distributed within a small 
area on a two-dimensional (2-D) disk. M t = 30 antennas 
transmit independent QPSK waveforms. The carrier frequency 
is 8.62 GHz. A maximum of L = 512 snapshots are consid- 
ered at the receive node. The received signal is corrupted by 
zero mean Gaussian noise. The SNR is set to 20 dB. There are 
two targets located at 6*i = — 3°and6>i = —2°, with reflection 
coefficients /3 k = 1, k = 1, 2. A jammer is located at 0° with 
an unknown Gaussian random waveform and with amplitude 
10, i.e., 20 dB above the target reflection coefficients f3 k - We 
sample the angle space by increments of 0.2° from —5° to 5°, 
i.e., a = [-5°, -4.8°, . . . , 4.8°, 5°]. We compare the perfor- 
mance of the proposed method and three approaches 12, i.e., 
the Capon, APES and GLRT techniques. 

Figs. [Ushows the moduli of the estimated reflection coef- 
ficients j3 k , as functions of the azimuthal angle for (a) M r = 1 
and (b) 10 receive antennas, respectively. In both (a) and (b), 
the top three curves correspond to the azimuthal estimates ob- 
tained via Capon, APES and GLRT, using 512 snapshots. The 
bottom curve is the result of the proposed approach, obtained 
using 15 snapshots only. One can see that in the case of us- 
ing only one receive node, the presence of the two targets is 



clearly evident via the proposed method based on 15 snap- 
shots only. The other methods produce spurious peaks away 
from the target locations. When the measurements of multi- 
ple receive nodes are used at a fusion center, the proposed ap- 
proach can yield similar performance to the other three meth- 
ods. However, the comparison methods would have to trans- 
mit to the fusion center 512 received samples each, while in 
the proposed approach, each node would need to transmit 15 
samples each. 

The threshold /j in © affects DOA estimation for the pro- 
posed method. The increase in /j, within a range can reduce 
the ripples of DOA estimates at the non-target azimuth an- 
gles at the expense of the accuracy of the target-reflection- 
coefficient estimates. The increase in [i can also reduce the 
complexity of ([8]l because the constraint is looser than that 
of smaller /j,. If /j, is too large, however, the £i-norm mini- 
mization does not work. In FigQ] a relatively large threshold, 
i.e., /i = 3, was used for the single receive node case. As a 
result, the CS method yielded less accurate estimates of the 
reflection coefficients magnitude than the Capon and GLRT, 
but with very few ripples. 

Fig. [2] shows the effect of the number of snapshots L 
on the DOA estimates of the proposed method and the other 
methods. Here, we consider the case of one receive antenna 
only. In order to quantify the performances of DOA estima- 
tion, we define the ratio of the square amplitude of the DOA 
estimate at the target azimuth angle to the sum of the square 
amplitude of DOA estimates at other angles as the peak-to- 
ripple ratio (PRR). As shown in Fig. [2], although the increase 
in L can improve the PPR of these four methods, the increase 
is much faster for the CS method. 

6. CONCLUSION 
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A compressive sensing method has been proposed to estimate 
the DOA of targets for MIMO radar systems. The DOA of tar- 
gets can construct a sparse vector in the angle space. There- 
fore, we can solve for this sparse vector by ^-norm mini- 
mization with many fewer samples than conventional meth- 
ods, i.e. the Capon, APES and GLRT techniques. The pro- 
posed method is superior to these conventional methods when 
one receive antenna is active. If multiple receive antennas are 
used, the proposed approach can yield similar performance to 
the other three methods, but by using far fewer samples. 
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Fig. 1. DOA estimation using (a) 1 receive antenna and (b) 10 receive antennas. In 
both (a) and (b), the top three curves were obtained using 512 snapshots. The bottom 
curve was obtained using 15 snapshots only. 
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Fig. 2. PRR using one receive antenna. 



